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Abstract. We use the classical umbral calculus to describe Riordan arrays. Here, a Riordan array is 
generated by a pair of umbrae, and this provides efficient proofs of several basic results of the theory such 
as the multiplication rule, the recursive properties, the fundamental theorem and the connection with 
Sheffer sequences. In particular, we show that the fundamental theorem turns out to be a reformulation of 
the umbral Abel identity. As an application, we give an elementary approach to the problem of extending 
integer powers of Riordan arrays to complex powers in such a way that additivity of the exponents is 
preserved. Also, ordinary Riordan arrays are studied within the classical umbral perspective and some 
combinatorial identities are discussed regarding Catalan numbers, Fibonacci numbers and Chebyshev 
polynomials. 



1. Introduction 

An exponential Riordan array is an infinite lower triangular matrix generated by two formal exponential 
series A{z) = J2n>o^" fir ^^'^ ^(^) ^ Sn>i ^" frf' ^'^^^ oo 7^ and bi ^ 0, such that, denoting this 
matrix by {A, B), its (n, fc)-entry is given by 



Riordan arrays were introduced by Shapiro et.al. in i26: as a generalization of the well known Pascal 
triangle and other types of combinatorial arrays (Rogers also studied in |20) equivalent generalizations of 
the Pascal triangle under the name of renewal arrays). Note that the generators A and B used in [26] are 
ordinary formal power series instead of series of exponential type. In this case the {n, fc)-entry of {A, B) is 
given by {A, B)n.k — [z"] {A{z) _B(z)'')). Riordan arrays form a group under matrix multiplication. They 
have been extensively studied and characterized in connection with combinatorial identities, recursion 
properties and walk problems (see for instance [5l [Ml [25l [26l I28j l. 

The purpose of this paper is to give a promising symbolic treatment of the exponential Riordan group 
based on a renewed approach to umbral calculus initiated in 1994 by Rota and Taylor [23]. In the last 
decade, this approach has been continued mostly by Di Nardo, Niederhausen and Senato [H [TOl [11] and 
more recently by PetruUo [9l I17[ I19j . We will refer to this new version as the classical umbral calculus, to 
distinguish it from the more established treatment of umbral calculus using operator theory, as presented 
for instance in Roman's book [21) . Under the old point of view of umbral calculus, Riordan arrays are 
also known as recursive matrices (see Barnabei, Brini and Nicoletti [2])- 
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The classical umbral calculus consists of formal operations over symbols, called umbrae and usually 
denoted by greek letters, whose powers represent sequences of numbers. A main feature of the classical 
umbral calculus is that the same sequence of numbers can be represented by distinct umbrae. That allows 
to suitable encode a broad family of generating functions. Nonetheless, more than an elegant method 
to encode and describe well-known concepts involving generating functions, the renewed umbral setting 
provides remarkable computational simplifications and conceptual clarifications in several contexts. For 
instance, let (s„(x)) be a Sheffer polynomial sequence, classically defined by (see [21]) 

1 + Vs„(a;)— =A(z)e-^W. 
^ — ' n\ 

n>l 

A key observation is that a sequence (s„(a;)) is of Sheffer type if and only if 

n 

(1) S„(X) =^(AB)n,fc^^ 

fc=0 

Formula ([T]) shows that Riordan arrays essentially provide a different point of view on Sheffer sequences. 
In [7], Sheffer sequences are described in terms of umbrae and, as an application, an umbral represen- 
tation for Riordan arrays is derived in [S]. In this work we elaborate more on Riordan arrays using a 
normalized version of such umbral representation (see Definition [TS] in Section 13.11 and compare with 
Corollary 5.5 in 8J. Tipically, the type of Sheffer sequences used in applications are those associated to 
classical orthogonal polynomials, where so{x) — 1. Taking this fact into consideration, we concentrate on 
normalized Riordan arrays; that is, those corresponding to monic Sheffer sequences (so(x) = 1). This is 
not a great loss of generality, as lots of the arrays considered in enumerative problems (Pascal, Stirling of 
first and second kind, and many others) satisfy the condition [A, B)n^n — 1; that is ao = 1 = 6i. Besides, 
often the Riordan group is defined for normalized Riordan arrays (see for instance jT2l [26]). 

In the context of combinatorial identities, the classical Abel identity provides a deep generalization of 
the well-known binomial identity 6 . Its umbral counterpart is written as (see Section 12. 2|) 

(2) (7 + arc,J2 (fj (7 + fc.a)"-^-a(a - k.a)"-' . 

k=Q ^ ^ 

It turns out that the fundamental theorem of Riordan arrays (Theorem 13. ip is simply a reformulation of 
the umbral Abel identity. In particular, formula ([T]) is stated as (Theorem 13. 3p 

(3) (7 + x.l3.a^r ^ E (fc) + ^•")"~' ■ 

As an application of the fundamental theorem, we give a short umbral proof of the multiplication rule for 
Riordan arrays (Theorem 13. 2p . Note that all known results about Riordan arrays can be restated in the 
umbral language using generating functions. However, more than doing a plain translation, we emphasize 
that the goal of this paper is to derive these results working within the umbral syntax in the first place, 
underlying in this way the computational convenience of following the umbral approach. As an instance, 
we give an elementary treatment to the subject of generalized powers of Riordan arrays, which are used 
to describe one-parameter subgroups in combinatorial physics [3] [12] . 

This paper is organized as follows. In Section [2] we recall the main features of the classical umbral 
calculus and review the umbral Abel identity and the Lagrange inversion formula. In Section [3] we 
unfold our umbral approach to exponential Riordan arrays, putting special emphasis on the fundamental 
theorem, from which short umbral proofs can be given for the multiplication rule of Riordan arrays and 
for the isomorphism between the Sheffer and the Riordan groups. We also give an umbral description of 
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several important Riordan subgroups. The simplicity of the umbral notation allows us to give a short 
presentation of the main recursive properties of Riordan arrays in Sectional In Section [SJ a few Riordan 
arrays are shown for which we make explicit some of the features already discussed. In Section |6l we give 
an elementary approach to generalized powers of Riordan arrays that recovers some of the results in [12j . 
Section [7] is devoted to ordinary Riordan arrays. Finally, in Section [8] we apply the umbral syntax to 
the Catalan triangle introduced by Shapiro [24 . Moreover, we recover some known facts and point out 
connections with ballot numbers, Fibonacci numbers and the Chebyshev polynomials of the first kind. 

2. The classical umbral calculus 

This section briefly recalls the basics of umbral calculus underlying much of this paper. Umbral Abel 
polynomials and the umbral Abel identity are also reviewed. For more information on this we refer the 
reader to [lOl [TTl [H [19] . 

2.1. Preliminaries. The basic ingredients to perform classical umbral calculus are: 

(i) A commutative integral domain R with identity 1 (we will usually set i? = C or i? = C[a;, y] for 
simplicity). 

(ii) An alphabet A = {a,/?, 7, . . .}, whose elements are called umbrae. 

(iii) A linear functional E : R[A\ — >■ R called evaluation such that 

• E[l] = 1 and 

• E[x'^y"^ a^j3^ ■ ■ ■ 7'''] = a;"y™ i?[Q;*]i?[/?^] ■ • • -^[7*''] (this property is called uncorrelation) . 

(iv) Two special umbrae: e [augmentation) and v [unity) such that £^[e"] — (5o,n and ~ 1 for 
aU 71 > 0. 

An umbra a represents a sequence (a„) = (a„)„£iN in R in the sense that i?[Q;"] = a„ for all n > (this 
forces flo = 1, which is not a great loss of generality). We refer to £'[a"] as the n-th moment of a and 
we assume that each sequence (a„) is represented by infinite many distinct umbrae, sometimes denoted 
by a, a' , a", . . . (this property is named saturation). An element p e R[A\ is called umbral polynomial 
and we define the support of p to be the set of all umbrae that appear in p. Two umbral polynomials 
p and q in R[A\ are said to be uncorrelated if E[pq] = E[p]E[q]. Therefore, distinct umbrae are always 
uncorrclated. On A we define two equivalence relations: 

Umbral equivalence: a ~ 7 if and only if E[a\ = -E[7], and 
Similarity: a = 7 if and only if _E[a"] = -^[7"], for all n > 0. 

Both ~ and = are extended to the whole of R[A\ in the obvious way. A very useful device in the classical 
umbral calculus is the generating function of an umbra. Given an umbra a we may consider the formal 
exponential series 

00 „ 

(4) e"" := 1 + y a" ^ G R[A] [[z\] . 

n=l 

The evaluation E can be linearly extended to such a formal series so that 

(5) Uz) := E[e"^-] + ^ ^ ■ 

We shall write e"^ ~ /q(-z) to mean E[e°'^] = fa{z), and we call fa[z) the generating function (g.f. for 
short) of a. The g.f. of an umbral polynomial p is defined in an analogous way. Also, note that a = 7 if 
and only if e"^ ~ e'*'^. The table below gives the moments and the generating function of some important 
umbrae that will be useful throughout the paper. 
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name 


umbra 


g.f. 


moments 


augmentation 


e 


1 


1,0,0,... 


Bernoulli 


L 


z 

- 1 


1, 6i, 62, ■ • • (Bernoulli numbers) 


unity 


V 




1,1,1,... 


singleton 


X 


1+ z 


1,1,0,... 


Bell 


P 




1, i?2, i^a, . . . (Bell numbers) 


boolean unity 


V 


1 


1,2!,3!,... 



A key feature of the renewed approach to umbral calculus is given by the dot product. By starting 
from elements in the dot product provides new symbols called auxiliary umbrae. Henceforth, when 

a new family F of auxiliary umbrae is defined, we will extend the alphabet A so that F (Z A. Often, the 
umbrae of A are used to make explicit the action of the evaluation on the powers of an auxiliary umbra. 
For instance, the action of E on the dot product 7. a is defined by 

(6) f^.a{z) = f^{\0gfa{z)) . 

It readily follows from ([6]) that j = r] implies 7. a = rj.a and a. 7 = a.rj, for all a. In particular, we have 
e.a = e = a.e and v. a = a = a.v, for all a. Also, /3 and x satisfy x-P = P-X = The dot product k.a, 
with fc e Z, has g.f. fk.a{z) = fa{z)^- If fc > and if a' ,a'\ . . . , a'*' are k pairwise distinct auxiliary 
umbrae defined to be similar to a then we have k.a = a' + a" + • ■ • + a**"' , so that k.a can be thought of 
as a generalization of the usual notion of multiple of an umbra. The umbra — l.a is named the inverse 
of a and satisfies a + (—l.a) = 0. The dot product a.x is called the a-factorial umbra (see [llj). Since 
/a.x(^) = /a(log(l + z)) ^ e"'°g(i+^) = (1 + z)" = 1 + E„>i(«)nS' thcu we deduce 

(7) (a.x)"^(a)„, n>l, 

where («)„ = a{a — I) ■ ■ ■ (a — n + I). Analogously, since v = — l-x- — 1 and a.c = ca, it follows from 
that (q.w)" ~ (a. - l-x- - 1)" ^ (-l)"(-a)„, and finally 

(8) (a.u)" ~ (a + n-l)„, n > 1. 

In order to manipulate the formal composition of two g.f.'s we iterate the dot product by obtaining 
auxiliary umbrae of type 77.7.0;. Hence, the g.f. of 7./?. a is 

f-r.p.aiz) = A(l0g/^(l0g/„(z))) = ^(/a(z) - l) . 

We write a*"^* to refer to a given symbol in A, if it exists, satisfying a^^'-Kfi.a = x = a.(3.a^^'-K We call 
a*^^* the compositional inverse of a. Note that E[a] ^ for a*"^* to exist. The table below translates 
some operations among umbrae in terms of their corresponding g.f.'s. 



umbrae 


g.f.'s 


a + 7 


Cauchy product 


a 77 


Hadamard product 


0.7 


log composition 


a./3.7 


composition 
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2.2. Abel polynomials and the Lagrange inversion formula. Let a be an umbra. The derivative 
of a is an auxiliary umbra whose powers satisfy ~ rta''^^ for n > 1. In particular — 1 

and this implies that has compositional inverse ol'-^^'' . The derivative umbra provides a binomial-like 
expansion of the powers (7./3.a^)" in terms of the umbrae a and 7. In fact, we have {fa^{z) ~ 1)" = 
z"faiz)"- ~ 2"e("-")^ and then 

e(..^.«.)^ ^ - E ^'^ ^ ^'"-"^^ - E E^" ^Ir • 

m— m— fc— 

Therefore, by comparing the coefficients of |j we recover 

(9) (7./3.aJ"^E(^)7'=(fc.ar-'=. 

fe=0 ^ ^ 

The umbral Abel polynomials a„(7, a) = 7(7 — 71.0)""^ are straightforward generalizations of the classical 
Abel polynomials x{x — na)"~^. They are powerful tools that allows one to understand the connections 
between an umbra and its compositional inverse. Such polynomials were first studied (in the special 
case 7 = x) by Rota, Shen and Taylor [22: in connection with polynomials of binomial type. More 
recently, they have been applied to the cumulant theory by Di Nardo, Petrullo and Senato [3, and 
further generalized by Petrullo [171 UHl HI] . Abel polynomials satisfy the Abel Identity [6] , 



(x + vT = E (fc) iv + kaT-'^xix ~ ka) 
It generalizes in a direct way to umbral Abel polynomials as follows 



fc-i 



(10) (7 + a)" ^ E it) + k.ar-^a{a - k.a) 

One of the main consequences of (jlOp is a one-line proof of the Lagrange inversion formula. To this end, 
it is useful to introduce the umbra ^ which satisfies 

(11) J?!;„ ~7(7-n.a)"-\ n>l. 
Theorem 2.1 (Lagrange inversion formula). For all umbrae a, 7, we have 

(12) (7./3.a;-^>)"~ 7(7 -«•«)""', n>l, 

Proof. By applying then (ITT|) . and finally (ITUl) . we have R^^^.p.a^ = 7, that is K^^^, = J-P-a^^^K □ 
The translation of (|12p in terms of generating functions gives the familiar identity 

n[z"]/,((z/„(z))<->) = [z"~i]/;(z) 

\Ja\Z) ^ 

and the reader would be convinced that (I12[) is exactly the umbral coding for the Lagrange inversion 
formula. Thanks to equivalence we may restate the Abel identity PH)) in a different, but equivalent, 
way. Indeed, set rj = a.f3.a''^^K Hence, we have a = rj.fi.a^ and, via (HH), ?y" — <j{(t — n.a)"~^. This 
way, the Abel identity (flUl) can be rewritten as 

(13) (7 + ^./3.«J" ^ E it) + 

fc=0 ^ ^ 

By means of ^ and it is easily checked that 



(14) (a-|-7.^.a^)^ = 7j,./3.a^. 
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Now, let £^ „ denote an umbra satisfying £^ „ = —1.^^^, and then set :— M.^^ and := The 
umbrae ^„ and play a central role in the umbral theory of cumulants [Hllin]- Since £„ = — l.a./3.a^"^\ 
via ([T4)l we obtain x = £t> = (a + a^) = {2^)^.(3.a^. This says that 

(15) ^a;->, 

hence £f;^ = a, and this explains why is named the Lagrange involution of a. Such an umbra is also 
useful to recover a reformulation of ([T2|) that will reveal useful for our aims. Let fc > 1 and denote by (5**°' 
an umbra with moments (J**"')" ~ ^fc,n, for n > 1. By virtue of (ITSI) we have (5''''./3.(£„)^ = ^^'''./S.a^^^^ 
Furthermore, we use dH) and ([12]) to get (^) (fc.£„)"-'^' ~ (-n.a)"~'=, and then 

(16) n(fc.£„)"-'= ~ fc(-n.a)"-^ 

By taking into account ££^ = a we may exchange a and £„ and obtain 

(17) n{k.a)"-'' - ^(n.J^J"-^ 

Again, we stress that ((TB]) is the umbral coding of the following identity on formal power series, 

3. Exponential Riordan arrays 



In this section we use the syntax developed in Section [5] to give an umbral coding of some fundamental 
facts of the theory of the (exponential) Riordan arrays. Umbral expressions for several important Riordan 
subgroups are given while the group axioms are efficiently translated in terms of the umbral algebra. 
Henceforth, we set R — G[x,y]. 

3.1. The exponential Riordan group. 

Definition 18. Given two umbrae a and 7 we denote by (7,0;) the infinite lower triangular complex 
matrix whose entries {'y,a)n,k satisfy 

(19) (7, a)n,fe - Q (7 + k.ar-" for n, fc > . 

We name (7, a) the (exponential) Riordan array of a and 7. Riordan arrays (j>,q), for p,q E R[-A], are 
defined in an analogous way. 

Riordan arrays are often introduced as a generalization of the Pascal triangle. Through ([T9l) the 
umbral syntax highlights this circumstance even more. Riordan arrays act on a column array of complex 
numbers in the obvious way. We write (7, a) a = b to express that the column array b is obtained by 
multiplying the column array a by (7, a) on its left. Since we deal with a saturated umbral calculus, then 
there exist umbrae rj and a; whose moments are the entries of a and b respectively. Hence, we may think 
of the Riordan action as an action on the alphabet A, and write (7, a) • = w instead of (7, a) a — b. 

Theorem 3.1 (Fundamental Theorem). For all umbrae a,j,r] we have 

(20) (7,0;) -T? = 7 + ?7./3.a^. 
Proof. By virtue of (fT9|) we have 



((7, a) 'vT^Y. il) + k-ar^'v", n>0. 

1 — n V / 



fe=0 

and the claim follows via (fT3l). □ 
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Then, from an umbral point of view the Fundamental theorem of Riordan arrays is nothing but a refor- 
mulation of the Abel identity pUj) . One of the main consequence of this theorem is obtained by setting 
rj = 1. In fact, in this case the umbra 7 + p.a^ represents the column vector whose nth entry is the sum 
of all entries occurring in the nth row of (7, a), 

n 

(21) (7 + /3.a,)"^^(7,a)«,fe. 

fe=0 

Consider two Riordan arrays (7, a) and (cr, p). We denote by (7, a){a, p) the array obtained by multiplying 
(7, a) by (ct, p). Hence, the following theorem can be stated. 

Theorem 3.2. The set £Hio of all Riordan arrays is closed under matrix multiplication. In particular, 
we have 

(22) (7, a)(cr, p) = (7 + cr.^.a^ , a + p./S.a^) . 
Proof. From and (fH|) we have 

(7,Q;)(cr,p) •r/ = (7,0;) • {a + i]./3.p^) = (7 + cr./3.a^) + T].f3.{a + p./S.a.^)^ . 
We may conclude that (7, a){a, p) •rj= {-^ + a.p.a^ , a + p.fH.a^) • rj for all 77, and then ([^ holds. □ 

The array (e, e) is the unity with respect to the matrix multiplication. Moreover, thanks to Theorem 1 2. 1[ 
and by virtue of we may deduce that the inverse of (7, a) is given by 

(7,a)-i = (£,.„,£„). 

Hence, *Hio is a group with respect to matrix multiplication and we name it the (exponential) Riordan 
group. 

3.2. The Sheffer group. 

Definition 23. Given umbrae a and 7, the polynomial sequence (s„(x)) defined by 

(24) 5n{x) + x.jS.a^Y for n>0, 
will be named the Sheffer sequence of (7, a). 

Sheffer sequences are classical polynomial sequences. They include several classes of orthogonal systems: 
Hermite, Charlier, Laguerre, Meixner of the first kind, and Meixner of the second kind. Sheffer sequences 
have been deeply studied within finite operator calculus by Roman |21j . Applications of umbrae to Sheffer 
sequences are in p!51 [55] , and more extensively in 18^ . By comparing ([50]) and (IM)) the following relation 
among Sheffer sequences and Riordan arrays comes trivially. 

Theorem 3.3. The Sheffer sequence (s„(a;)) 0/(7,0;) is represented by (7,0;) • x. Equivalently, the 
following expansion is possible 

n 

(25) s„(a:) = ^ (7,a)„,fca;'' . 

k=0 

Hence, by virtue of (|25p the following fact is achieved: the Riordan array (7, a) is nothing but the array 
of coefficients of the Sheffer sequence (s„(a;)) in terms of the canonical basis (x"). Note that, from ([501) 
we obtain 



(26) 



(7, a) • (ry + cj) ^ (7, a) • rj + {e,a) • uj 
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If (s„(a;)) and (p„(a;)) are the ShefFer sequences of (7,0?) and (e,a) respectively, by setting rj — x and 
a; = zy in (j26p . and by taking the nth moment of both sides, we may recover the well-known Sheffer 
identity, 

k=0 ^ 

Sheffer sequences are endowed with a binary operation called umbral composition. More precisely, if 
(s„(x)) and (r„(a;)) are the Sheffer sequences of (7, a) and {a, p) respectively, then the umbral composition 
of (s„(a;)) with (r„(a;)) is the polynomial sequence (sr„(x)) defined by 

n 

(27) sr„ {x) = ^ (7, a)n,k U {x) for n > . 

fc=0 

By expanding the r„(a:)'s according to (|25l) we may check that the coefficient of x'^ in sr„(a;) is the 
(n, fc)-entry of the Riordan array (7, a){a, p). This way, from (|22l) we obtain that (sc„(a;)) is the ShefFer 
sequence of {-^ + a.(3.a^,a + p.(3.a^). Also, we deduce that (s„(a:)) and {xn{x)) are inverses of each other 
if and only if (cr, p) — (£.^, and that the Sheffer sequence (x") of (e,£) is the unity with respect to 
the umbral composition. We can conclude that the set ©f)eff of all Sheffer sequences is a group under 
the umbral composition, which we name the Sheffer group. The following theorem is now obvious. 

Theorem 3.4. The Sheffer group and the Riordan group are isomorphic. 

3.3. Riordan subgroups. We give umbral characterizations for some important Riordan subgroups. 

The Appell subgroup. We denote by Appell the set of all arrays of type (7,e). Given any two elements 
(7,e) and (cr, e) in Appell, it follows that 

(7,e)(cr,e) = (7 + cr,£) and (7, e)"^ = (-I.7, e) . 

This shows that Appell is a subgroup. We can also see from the left formula above that it is an abelian 
subgroup. Being = x and = v, then Appell acts on the alphabet of umbrae following 

(7, e) • r/ = 7 + r/ . 

Sheffer sequences relative to the Appell subgroup are named Appell sequences. 

The Associated subgroup. The collection of arrays of the form (£,a) is denoted by Assoc. Given any 
(e,a) and (e,p) in Assoc, we have 

{e,a){e, p) = {e,a + p.l3.a^) and {e,a)^^ — {e, 2^) . 

The action of Assoc on A is encoded by 

{e,a) •■q = Tj.p.a^ . 

The corresponding Sheffer sequences are named associated sequences or binomial sequences. 

The Bell subgroup. Of particular interest are arrays of the form (a, a). We denote the collection of all 
such arrays by Bell. Given any {a, a) and (cr, cr) we have 

(a, a)(cr, cr) = (a + CT./S.a^,, a + CT./3.a^) and (a, a)^"'" = (£„,£„) . 

The Be/Z-action is encoded by 

{a,a) • -q = a + r].(3.a^ . 
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The Stabilizer subgroups. Given any 77 £ ^, the stabilizer Stab{r]) of 77 (with respect to the ?lio-action) is 

Stab{rj) = I (7, a) £ Ulio : ^ + ri.p.a^ = Tyj . 
Since (7,0!) G Stab{ri) if and only if 7 = — l.r/./J.a^, then we have 
Stab{e) = I (7, a) e ^io : 7 = e| = Assoc. 

Stab{v) = I (7, a) G ^Hio : 7 = 1 — l./3.a^| = Stock (Stochastic subgroup). 

Stabix) = {(7,a) G$Hio : 7 = x-l-a,,}- 

We give explicit umbral representations for the {n, fc)-entry of Riordan arrays belonging to some of the 
subgroups mentioned above. The stabilizers case can be recover by applying ([T3l) and (fT9|). 



Subgroup 




Appeh (7,e) 




Associated (e, a) 


(:)(*.«)•-' 


Bell {a, a) 




Stab{ri) (rj ~ l.Ty./J.a^,, a) 





4. Recursive properties of Riordan arrays 
In this section the main recursive properties of a Riordan array (7, a) are established. 
Theorem 4.1. The Riordan array (7, a) satisfies the following recursive properties: 

(28) 

and 

(29) (7,a)„,fc 

Proof. From we obtain 



(7, a)„,fe - ^ X! j ^ ) (7: a)n-i-i,fc-i fc > 1 , 

i=0 ^ 



n -^—^ fk — 1 + i\ 

I J ) -He (77 ")«-!, fc-i+i n,k>l 
1=0 ^ 



(7, «)„,. ^ (^) (7 + (fc - l).a + a)-' Q g (7 + (fc - l).a)"-'=~'^ , 

and then ((28|) follows thanks to elementary computations. Moreover by applying the Abel identity ([10 
we have 

(7,a)",fc 



(7 + (fc - !)•« + a)" ~ Q II J ^) (7 + - 1 + 'a{a- i-af-^ , 



and (I29p is a consequence of (fTTj) . 



□ 
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Equation (j28p expresses the (n, fc)-entry of (7, a) as a linear combination of the entries in the preceding 
column up to the preceding row. Equation (|29p gives the (n, fc)-entry of (7, a) as a linear combination of 
the entries in the preceding row starting from the preceding column (note that from column n — 1 on all 
entries are zero). This result is the umbral equivalent version of Rogers' recursion formula for ordinary 
Riordan arrays [20] (see also [28]). The sequence of moments of is the A-sequence of (7, a). Since 
^a-l3-a^ = a and ~ a{a — n.o;)""^, a further recurrence relation can be derived from (|13p as follows. 

Theorem 4.2. The Riordan array (7, a) satisfies the following recursive properties: 

/ \ n — k 

(30) (7, a)n,k ^ i^j (^•^°)' (T' n,k>l. 



Proof. From p9|) and (|T3|) we obtain 



□ 



5. Examples 

In this section we show a few examples of classical arrays that are in CHio. We derive some of their 
recursive properties and make explicit the related ShefFer sequences. 

Identity array. The identity array is the Riordan array I — {e, e). Its Sheffer sequence is (x"). 

Pascal array. The Pascal array is the Riordan array P = {v, e). So we have Pn.k = {^)- 
(1) Note that P • v = 2 and then we recover the well-known identity 



2" = E 

fc=0 

(2) From P»a:; = l + a;we deduce that 

I — n \ / 



fe=0 

Stirling arrays. The Stirling array of the second kind is defined by S' = (e, — l.t). Its entries Sn,kS are 
named Stirling numbers of the second kind and are denoted by S{n,kys. Some fundamental properties 
of this array can be recovered as follows. 

(1) Since (— l.t)^ = v then from ([^(7| we obtain S • v = P and then the row-sum of this array gives 
the sequence of Bell numbers: 

n 
fc=0 

(2) Let ((/>„(a;)) be the binomial sequence represented by x.f5. The 0„(a;)'s are called exponential 
polynomials [21]. Since S • x = x.(3, we deduce that the Sheffer sequence of S is (0,1 (.x)), hence 

n 
fc=0 
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(3) From (^5)) we recover 



n — k 



Since ^-i.t = — l.t)^ ^\ and being (— l.t)^ ^> = v'- = x-Xj we obtain = — l.i-x- 

Equivalence ([^ gives 

fc) = ^ |]7^ ~ ^ + 5(n - 1 - fc - 1) , 
where the "^i's are sometimes called Cauchy numbers [TSl. Finally from ([50)1 we obtain 

with the '^^''''s denoting generalized Cauchy numbers, whose g.f. is the fcth power of the g.f. of 
the "Ti's. 

The Stirling array of the first kind is defined by s = S^^. Being = — l.^.i.^ = l-x, it immediately 
follows that s = (e, i-x)- 

(1) Via generating functions it can be checked that (t.x)x> = X-X ^^'^ then we obtain s»v = x- This 
way, as s(0, 0) = 1 = s(l, 1) and s{n, 0) = for n > 1 then we have 



— s(n, k) for n >2 . 

k=0 

(2) From s • x = x.x we deduce that 

n 

{x)n = ^ s(n, k) x" . 

fc=0 

and then the Sheffer sequence of s is the sequence of falling factorials ( (2^)11) 

(3) Note that (x-x)" ^ - 1)!. Then, from ^ we recover 

n — k 

' 1 + 



i=0 

Since ^^.^ = l then ((29)) gives 

n — k 

k 



^) - I E ^ (-1)' - 1 - fc - 1) . 
A; I + 1 



("'^) = rE \b,s{n-l-i,k-l). 



Finally from ([30]) we obtain 

y \ n — k 

with 6^*°' denoting a generalized Bernoulli number [B]. 
We conclude this example by noting that s • (x.x) = {s ■ • x = x. This says 

n 

k=a 

and then s stores the transition matrix from ((a;)„) to (x"). 
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6. Generalized powers in the Riordan group 

Let j E K and (7, a) £ Uiio. In this section we give an elementary approach to the problem of extending 
the power j in (7, a)^ to any complex number c, such that (7, a)^ £ Ulio. This subject has been already 
discussed in [3| I12j. where some physical insights are underlined. In our context, by (7,0)^ we mean 

(31) i-f,ay = (7<'',a<'')---(7*'^""')' 

" V ' 

i 

where the 7''' and a''' are respectively umbrae similar to 7 and a, for all 1 < i < j- By denoting 
(P.a^)*^ = ^.a^'.^.a^-^'. • • • ./3.a^'' and setting {l3.a^)*° = v, we can write dSj) as 



(32) (7,a)^- = K]7<"-(/3.aJ*(^-^\£"*"-(/5-«- 



(t-i) 

i=l / 



Formula p2p makes explicit the two umbrae this is a Riordan array of, so that we can write (7, a)^ G 9iio. 
In general, given an arbitrary c S C, in order to say that (7,a)'^ G *Hio, we need to find umbrae a and 
p such that (7,a)^ = (cr, p). To achieve this goal, denote by (7, a) + (c, p) the componentwise matrix 
summation. Note that (7,0) + (a, p) ^ (7 + <t, a + p) and in general (7,0) + (u, p) is not an element of *Kio . 
For example, the matrix represented by (7, a) — (e, e) is not in dKxo as the (n, fc)-entry ((7, a) — (e, e))^ j, 
equals if n = fc, and reduces to (7, Q;)„^fe otherwise (note that (7,0;) G fHio implies that (7,q;)„^„ ~ 1 
for all 71 G M). Nevertheless, let us set 

(33) (7,a)-:=^(')((7,a)-(e,e))'. 

It readily follows that by choosing c = m G W the generalized power f33)) matches the usual power 
(7, a)™. However, it is not immediately clear that (7, a)^"^ G $Hio for a generic c G C (it is not even clear 
that (7, a)^^ is well defined). The next step is to find two auxiliary umbrae 7^^^ and a^'^ so that 

(7,a)- = (7^a*=) 

Obviously, if such umbrae 7^'^ and a^'^ exist, they must be given in terms of 7 and a. Formula p2p 
suggests that the umbra 7^^^ should depend on 7 and a, whereas the umbra a^'^ should depend only on 
a. Now, for j > 1 we have 

((7, a) - (£,£))' = ((7"',««) - •• • ((7<^a'^') - , 

3 

from where we get the following expansion. 



fc<ii <---<ij_i <n 

By setting Iq = k and = n and applying (|19l) . we obtain 

((7,«)-(e,e))' E f, , " / , V-T'^'+^^-""')'^"^-^---(V^'+^o.a<^')'^"": 



where we use the multinomial coefficient 



Ij — Ij-i, . . . ,li — IqJ {Ij — Ij-iY- ■ • • (^1 — ^0)! 
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Therefore, from (I33p we can write 

]>1 ^-^^ la<li<---<lj-i<h ^ ^ 

Note that for j > n — k, the sum above is zero. This shows that the entries (7, a)^'^^. involve finite 
summations for any n and fc; hence (7, a)*'^ is well defined. In particular, when n = k we get (7, a)*'^„ ~ 1 
for all n > 0. Setting k = 0, we define the auxiliary umbra 7^'^ as the umbra having moments 

(34) (7^T^(£,£ko + E(-) E L-/.'' )(7"'+^,-i.a«')"-'-^---(7'^>y^- 

Analogously, setting k — 1 and taking 7 = e, we define the auxiliary umbra a^'^ as the umbra such that 
for n > 1, it satisfies 

(35) (a*T-^^(^,^)n,i+Ef') E f, , " , ,1 
We emphasize that 7''''^ and a^"^ are well defined umbrae. Therefore (7^^^, a^'^) S $Hio and by (|19p. we have 

(y^aH..-Q(y^ + fc.«*^)""^ 

Note that (7,a)^';o — (7'''^)" ~ {^^'^ , a^'^)n.o- In fact, more is true. 

Theorem 6.1. For any c G C and (7, a) £ U\io, the generalized power (7, a)^"^ defined in (j33p is m *Hio. 
Moreover, 

(7,a)- = (7^«*^), 
where 7^*^ and a-l- are given by p4p anrf (I35p . 

Proof. Let j G IN. Using recursively Theorem 13.11 in particular Formula ()13p . we verify that 

E7'"-(/3-a^)*^'~'M - (7^')" and [E - fo^' ^ > 0. 

Therefore, for all j G IN we have (7, a)*-' = (7, «)■' = (7^"' , a-'--' ) . Moreover, observe that {'y,ct)nk ^^"^ 
(f^)E[{'-f^'^ + k.a^'^)^^''] are polynomials in c that agree for all c — j G M. Hence, they are equal. □ 

In particular, when a = e we have 7^^^ = c.7 and e'^'^ = e, so that we can write (7,6)*'^ = (c.7,e). 
Furthermore, for any ci, C2 G C note that 

(7,e)^^'''^''''' = ((ci +C2).7,e) = (ci.7 + C2.7, e) = (ci.7, e) (02.7, e) = (7, e)*''' (7, e)^''" . 

Indeed, the property above is true in general, as we can readily check that 

(7,ar^-^^=E^t'')((^'")-(^'^))' 



3>0 



^j>0 ^-^ ' / \j>0 

(7, a) (7, a) 
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7. Ordinary Riordan arrays 

In this section we present the umbral syntax to handle ordinary Riordan arrays. We name ordinary 
Riordan array an infinite lower triangular matrix [7, a] defined by 

(36) [^,«]„^,^____. 

Since n![7,a]„^fc = kl{-f,a)n,k then all results of Section [3] and the recursions in Section |4] can be easily 
adapted to these arrays. For instance, from the identity (|13|) we have 



and this explains how the Fundamental Theorem 13.11 applies to ordinary arrays. More explicitly, if 
[7, a]a ~ b then the kth entry of the column vector a is -^^f^, and the nth entry of b is ^^'''^"'"''^f ""^"^ ^ - 
Moreover, elementary computations shows that 

(7, a) = {a, p){t, tt) if and only if [7, a] = [a, p] [t, tt], 

so that the set of all ordinary arrays is a group under the usual matrix multiplication. The identity is 
[e, e] and the inverse of an array is determined by means of 

[7, = £J. 

If 1 has all the entries equal to 1 then [7, a]l stores the row-sums of [7, a]. Then, from -u" ~ n! we obtain 
the following analogue of (|2T|) . 

(7 + u./3.a^)" 

-2^[7,a]n,fc- 

k=0 

The rules ([SS)) . (f29| and (|30|) translate into the recursions below, 

n—k j 

(38) [7,a]„,/c ^ ^ — [7,a]n-i-j,fc~i n,k>l, 

1=0 

(39) [7, a]n.k - -f a]n-i,k-i+i n,k>l, 
and 

(40) ['y,a]n,k^22 — '—^h^Oi]n-k,i n,k>l. 

i=o 

As in the exponential setting we may define the ordinary Sheffer sequence of [7, a] to be the polynomial 
sequence (s„(a;)) such that 

n 

Snix) = y^J'^,a]n,kx''. 
k=0 

By applying ([37]) we obtain 

Sn{x) ~ n! (7 + v.x.p.a^)^ . 

The umbral composition of the Sheffer sequences of [7, a] with the Sheffer sequence of [a, p] is the Sheffer 
sequence of [7, a] [cr, p]. The identity sequence (x") is the Sheffer sequence of [e,e]. In particular, the 
ordinary Riordan group IHio' and the ordinary Sheffer group 6ficff are isomorphic, and both of them 
are isomorphic to their exponential analogue. Note that the isomorphism between TOo' and ©f)eff' is one 
of the main results in [13] (see also [30]). The ordinary Appell sequences and the ordinary associated 
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sequences correspond to arrays of type [7,2] and [s,a] respectively. Finally, the ordinary ShefFer identity 
says that 

n 
k=Q 

whenever z'^ = J2^^o x^u'^^^ for all fc = 1, 2, . . . , n, and with (p„(a;)) denoting the sequence of [e, a]. 

By applying ([8]) it is easily seen that the Pascal array can be also realized as the ordinary array [w, "O]. 
In the same fashion as in the exponential setting, where (pS)) . and ([501 reduce to trivial identities, 
the equivalences ((38)) . ((39)) and ((40)) provides well known recursive properties of the binomial coefficients. 
In fact, by comparing ([7]) and ((Sj) we have n{k.v)^^^ ~ k{n.x)''^^^ . By virtue of ((T7)) this provides = x 
and the recursive identities are 




Moreover, from £^ = — l.^ we obtain 

(44) = [-l.X,-l.x]- 

We now consider an important ordinary array that involves Catalan numbers. Let be an umbra with 
moments ^ ~ n!C„, n > 1, with C„ = denoting the nth Catalan number. From ()16|) we have 

(n + l)(£x.-i)" — ((" + l)-"*^)" — (2'^)n- This way, we have £x--i = then £^ = x- — 1 and finally 
= V. Via ((28)) . ((29)) and ((30| we obtain the recursive properties of the array [1?,^], 

n—k 

(45) [<r,?]n,fc = ^ C'^';, <;]„-i-i,fc-i, 

i=0 



n — k 
i=0 



K, ^ ?^-l,fc-l- 



Also, we have 



(46) [^,^]-i = [x.-l,X.-l]. 

Moreover, via ()17p and (IMl) we obtain an explicit formula for <^]n,k- In fact, from (n + l) ((fc+l).<;) 
(fc + l)((n + ~ (fc + l)(2n - fc)„ we have 

,^ , , , k + 1 f2n — k 

(47) [^,^]„,fc 



n — k 



n-\-l\ n 



16 



J. AGAPITO, A. MESTRE, P. PETRULLO, AND M.M. TORRES 



These numbers are often called the ballot numbers and have interesting combinatorial insights [T]. The 
first rows of the array are 



/ 1 
1 
2 
5 



14 14 9 



V 



Since w = -u^, by means of ([TT|) and ([H)) we have (<; + v.p.^^)^ = v.j3.^^ = ^^.p.q^ = Hence, the 
row-sums of <;] are given by 



k=0 

8. A Catalan triangle, the Fibonacci numbers and the Chebyshev polynomials 



Applications of Riordan arrays mainly concern combinatorial problems [4] [51 1251 128) . In this section we 
apply the umbral syntax to discuss some properties of the Riordan array [2.<^, 2.<j] whose combinatorial 
insights have been discussed in several papers [3 [MJ [57] . We will recover most of their main properties 
by underlying some connections with ballot numbers, Fibonacci numbers and the Chebyshev polynomials 
of the second kind. First of all, note that from (jTT)) we have <;]„+i_i — C„+i and that by virtue of 
with k = 1 we recover the well-known convolution formula for the Catalan numbers, 

n 

Cn+l — ^ ^ CiCn-i- 

In particular, this says that 2.^ = <; -I- v./3.<;^ so that 

[2.<;,2.?] = [?,^][C,i;]. 

Then, [2.<;,2.<j] can be obtained by multiplying [<;,<;] by the Pascal array [v,v] on its right, 

i 
k 



[2.,,2.,]„,fc = ^ [,,^] 



In particular, via (HH) and we obtain a nice expression for [2.<;, 2.<;] ^, 

[2.<j,2.^]-i = [-1.X, -1.X][X- - 1,X- - !]• 

However, a more explicit expression of [2.<;, 2.<j]^^ may be obtained as follows. Note that (2n.x)"^^ 2± 
n(2.<;)"~-^. Hence, by virtue of (|T7)) we obtain ^2.^ = 2.x and then £2.^ = —2.x. Now, by applying (|T7)) 
again, we have (n -I- -f- 1).'?)" ~ (fc -f l)((n + l).2.x)" and the entries of the array are 

r„ „ -, A; + 1 /2n + 2 
[2.c,2.^]„,fe = — - 

n + i \ n — K 

This allows us to compute the first rows of the array 



[2.^,2.^]- 



/ 


1 










2 


1 








5 


4 


1 






14 


14 


6 1 






42 


48 


27 8 1 




V 








■ / 
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Thanks to (|39p . ([55)1 and (^0]) we obtain the foUowing recursions, 

[2..;, 2.^]n,k = [2.^, 2.^]n-i,k-i + 2 [2.?, 2.<;]„_i^fc + [2.?, 2.<;]„_i,fe+i, 



i=0 



and 



2k 



[2.<j,2.^]„,fe = ^ I . j[2.<j,2.^]„_fe,,. 

i=0 \ * / 

Now, we are interested in the Sheffer sequence (s„(a;)) of [2.<j, 2.<;]^^ = [—2.x, —2.x] satisfying 



(48) 



[2.^,2.?] 



Si{x) 

\ ■■ J \ J 



X 
^2 



Since {-2{k + l).^)""'' ^ + k+ l)„_fe, we obtain 

(49) 



fc=0 



n ~ k 



Such a polynomial sequence plays an important role in mathematics. In fact, from the umbral equivalence 
n\Sn{x) ~ (—2.x + v.x.l3.{—2.x)j,y^ we have the g.f. 



i+$:.„(.)z"=-^ 



{l + z^-xz l-(a;-2)z + z2 

It is now easy to recognize that s„(x) — C/„ (^^^) with (C^n(2:)) denoting the sequence of Chebyshev 
polynomials of the second kind [6]. It is not difficult to see that (s„(x)) satisfies the three-term recursion 

Sn{x) = {X - 2)Sn-l{x) -Sn-2{x), n > 2, 

with initial values So(x) — 1 and Si(a;) ~ x — 2. Such a polynomial sequence specializes to integers 
whenever x is integer, and then (|48p specializes to interesting matrix identities. Since the change of 
variable {x — 2) i— (2 — x) is such that Sn{x) i-^ (— I)"s„(a;), then it is sufficient to study the case 
x-2>0. 

Example I (Periodic sequences). If a; = 2 then the recursion s„(l) = — s„_2(l), with initial values 
Bo(f) = 1 and Si(l) — 0, is solved by the period 4 sequence 1,0, —1,0. More explicitly we have 



/ 1 

2 1 

5 4 1 

14 14 6 1 

42 48 27 8 1 



\ / 1 \ 



-1 



1 



/ 1 \ 

2 

22 
23 
24 



J \ ■■ J \ ■■ J 
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When X = 3 we have s„(3) = s„_i(3) — s„_2(3), with So (3) = Si(3) = 1. It is not difBcuh to see that the 
recursion is satisfied by the period 6 sequence 1, 1, 0, —1, —1, 0. We have 



/I \ / 1 \ 

2 1 

5 4 1 

14 14 6 1 

42 48 27 8 1 



1 



-1 
-1 



These two cases are better understood by recalhng that [/„(cos?/) 



/ 1 \ 

3 
32 
33 
34 

V ■■ J 



sin{n4-l)y 

sin y 



(see for instance [6]). 



Hence, we have s„(2) = C/„ (cos ^) and s„(3) = Un (cos -1) and this explains the periodic sequences. 

Example 2 (Natural numbers). When a; = 4 we have s„(4) = 2s„_i(4) — s„_2, with So(4) ~ 1 and 
Si(4) = 2. We obtain s„(4) — n + 1 and then 



/ 1 

2 1 

5 4 1 

14 14 6 1 

42 48 27 8 1 



\ / 1 \ 

2 
3 
4 
5 



V 



V J 



( 1 \ 

4 
42 

43 

44 

V ■■ ) 



In other terms, we have 



(fc + If (2n + 2 



k=0 



n+1 



- k 



Such an identity was proved by several authors and it is the motivation of the combinatorial approach of 
Chen et al. [5]- 

Example 3 (Fibonacci numbers). When x — b the recurrence s„(5) — 3s„_i(5) — Sn-2, with So(5) = 1 
and Si(5) = 3, is solved by the subsequence of Fibonacci numbers s„(5) = F2„+2- Then, the following 
identity relates Catalan numbers and Fibonacci numbers, 



/ 1 

2 1 

5 4 1 

14 14 6 1 

42 48 27 8 1 

V ; 



\ / 1 \ 

3 

8 

21 
55 

V ■■ ) 



( 1 \ 

5 

52 

53 
54 

V ■■ ) 



Equivalently, 
(50) 



^k + l (2n + 2\^ 



fe=0 



By taking into account (|49p . we may invert (15(11) and obtain 

n 

F2„+2 = (-!)" E 



n + k + 2 
n — k 



n + k + 1 
n — k — I 



(-5)^ 
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Moreover, being F2n+3 = ^2(n+i)+2 ^ -p2n+2, then we deduce an analogous formula for the Fibonacci 
numbers of an odd index, 

-^~«H-i."|:((":!r)-(„!:!,)).-3)' 

Formulae relating Fibonacci numbers and binomial coefficients are already known (see for instance [4], 
where a reformulation of an identity of Andrews that is analogous to the formulae above is proved by 
means of Riordan arrays) . 
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